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Abstract 



Radial tidal forces can be compressive instead of disruptive, a possibility that is frequently 
overlooked in high level physics courses. For example, radial tidal compression can emerge in 
extended stellar systems containing a smaller stellar cluster. For particular conditions the tidal 
field produced by this extended mass distribution can exert on the cluster it contains compressive 
effects instead of the common disruptive forces. This interesting aspect of gravity can be derived 
from standard relations given in many textbooks and introductory courses in astronomy and can 
serve as an opportunity to look closer at some aspects of gravitational physics, stellar dynamics, and 
differential geometry. The existence of compressive tides at the center of huge stellar systems might 
suggest new evolutionary scenarios for the formation of stars and primordial galactic formation 
processes. 
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I. INTRODUCTION 



Tides are the manifestation of a gradient of the gravitational force field induced by a mass 
above an extended body or a system of particles. In the solar system tidal perturbations 
act on compact bodies such as planets, moons, and comets.- On larger scales than the solar 
system, as in a galactic or cosmo logical context, we frequently deal with tidal deformations 
(and eventually even disruptions) of a stellar cluster by a galaxy or in galactic encounters. 
The study of tidal fields in the most simple cases can be done by analytically approximating 
the potential field of a nonspherical mass distribution by a spherical (that is, a point mass) 
potential. This approximation sometimes leads to gross estimates and is inappropriate for 
the analysis of tidal fields within mass distributions. 

Stars can be considered as particles with a long mean free path, and the clusters they 
form as collisionless fluid dynamical systems. On galactic or cosmological scales a huge mass 
distribution can interact gravitationally with a smaller system that travels through it or is 
inside of it, without any direct contact taking place between its constituents. For example, 
think of a globular cluster passing through a disk or near the center of a galaxy, or a dust 
cloud orbiting inside a star cluster. The tidal field produced by the larger system, acting 
on the smaller system it contains, needs to be computed inside the mass distribution. Tidal 
forces can in these cases behave very differently than in the cases encountered in the solar 
system. 

Although numerical simulations are usually necessary to study an extended mass distri- 
bution, the analytic study of some representative, not too approximate models is sometimes 
possible. The effect of tides due to an extended mass distribution has been an active field 
of research in the last two decades, especially since high speed numerical computers became 
available. Typical cases of interest are galactic encounters,-^ globular clusters under the in- 
fluence of the galactic mass distribution,-!^ and Oort cloud perturbations by the galactic 
field." However, most of these investigations were on disruptive tidal forces. 

We are accustomed to think of gravitation as a force field that decreases rapidly with 
the distance from the field source. This dependence is not necessarily true inside a mass 
distribution producing the field where the force can increase with the distance from the 
center of mass. A trivial example is a spherical homogeneous mass distribution with radius 
R and constant density po- To an outside observer at distance r > R the force field is 
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F{r) = -GM/r'^, with M = po|7^i?^ and the tidal force dF{r)/dr = 2GM/r^ > 0, which 
because it is positive is responsible for the typical tidal bulges, that is, disruptive effects. 
If we move inside the same distribution, we obtain F{r) = —GM{r)/r'^ = —^nGpor, but 
dF/dr = — IvrGpo < 0. The tidal field has become compressive because inside the body the 
gravitational strength increases with the distance from the center of the sphere. 

A constant density fluid system in gravitational equilibrium is an unrealistic case. The 
effect of the compressive nature of tidal fields has been studied in the context of disk shock- 
ing of star clusters.-^ Fast encounters with the galactic disk or bulge cause gravitational 
shocks where the entire cluster expands, but with a short period of contraction immedi- 
ately after the shock. Another attempt to understand this effect was made by Valluri^ who 
studied compressive tidal heating effects. The effects of the mean tidal field of a cluster of 
galaxies on the internal dynamics of a disk galaxy traveling through it was discussed, and 
the disk experiences a compressive tidal field within the core of the cluster. Das and Jog^ 
investigated the same effect for compressed molecular cloud dynamics. They suggested that 
the compressive tidal field in the center of fiat core early type galaxies and ultraluminous 
galaxies compresses molecular clouds producing the dense gas observed in the center of these 
galaxies. 

The nature of the radial compressive tide remains largely unexplored. It is important 
not only for its physical effects, but because it is a fundamental aspect of gravitation itself. 
In the following we formulate a tidal field theory in which this effect arises naturally and 
investigate its relevance in some real systems such as stellar globular clusters and galactic 
centers. 



II. THE CONSTRUCTION AND MEANING OF THE NEWTONIAN TIDAL 
TENSOR 



The Newtonian tidal tensor is a useful tool for the study of tidal fields. The use of the 
tidal tensor for this purpose is not new, although not common. It has been employed for 
example in the study of the angular momentum growth in protogalaxies and more recently 
for investigatingthe origin of angular momentum in galaxies through tidal torquing (see for 



example, Refs. 



and 



lOl ). Here we analyze the tidal field inside a mass distribution and 



show how the latter naturally contains regions of radial tidal forces sign reversal. 
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Because the Newtonian tidal tensor is frequently employed in a somewhat obscure way 
in the literature, it is useful to recall some of its properties and see how it naturally leads 
to the common tidal textbook formalism. 

Once a coordinate system and a frame of reference is defined (typically with its origin 
at the mass distribution's center of mass) and given the mass distribution's gravitational 
force field F(R) (with the vector R vector at a particle of the tidally perturbed system), 
the tidal forces arise as the manifestation of the gradient VF(R). The tidal force induced 
by a mass distribution on another body or another system of particles is defined as the 
difference between the gravitational force that the mass exerts at that point, and the mean 
gravitational force (F) it exerts on the whole body or system 

Ft(R)=F(R)-(F). (1) 

It can be shown that the mean gravitational force acting on the tidally perturbed system 
is equivalent to the force acting at its center of mass. Rem- Hence without any loss of 
generality, we can write the tidal fields as 

Ft(R) = AF = F(R) - F(Re^). (2) 

Once the gravitational force field is expressed in terms of the potential $(R) as F(R) = 
— V$(R), Ft is completely determined if we introduce the second rank symmetric Newtonian 
tidal force tensor— 



52$ 



(3) 



with i,j = (1,2,3) and Xi = x, X2 = y, X3 = z; t^j is symmetric due to the conservative 
character of gravitational fields. If we use Einstein's summation convention, the components 
of the tidal force Ft in its differential form is 

dF 9^$ 

= = -d^/""' = 

The tensor r is the Jacobian of the gravitational force field F and the negative Hessian 
matrix of the scalar potential function $. These analytic properties of r give information 
about the maximum, minimum, and saddle points of the potential surface ^{x,y,z).— The 
negative of the trace of r gives Poisson's equation, which relates the potential $ to the mass 
density p 
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Poisson's equation is nothing more than a trace invariant quantity of the tidal tensor matrix. 

It is possible to formulate the gravitational potential $ through Poisson's equation, 
and therefore the force field of heterogeneous spherical, spheroidal, and disk-like mass 
distributions.^^ This formulation is not straightforward, especially for irregular and com- 
plex mass distributions. From linear algebra we know that the trace of a Hessian r of a 
local surface patch is an invariant quantity under coordinate transformations and that it 
equals the sum of its eigenvalues, Aj, that is, Tr(r) = Differential geometry defines 

twice the value of the mean local curvature (the mean value of curvature along the axes) 
as H = iTr(r), and the Gaussian curvature is defined by the product of the eigenvalues: 
K = Yl^Xi, or by the determinant det(r). Hence, there seems to be a relation between 
the Newtonian tidal tensor, the curvature of the potential surface it represents, and mass 
density and tidal forces, which is reminiscent of the connection between spacetime curvature 
and matter in general relativity. This relation is not a coincidence. 

The Newtonian tidal tensor Tij is the classical counterpart of the more general fourth 
rank Riemann spacetime curvature tensor of general relativity, which is why it is called 
Newtonian. Those familiar with general relativity know that the Riemann tensor, R^ji, is 
also known as the tidal force tensor. It can be showni^ that in the weak field limit the 
Riemann tensor is the classical Newtonian tidal force tensor, that is, -R'ooi = '^ij- Curvature 
in general relativity is not associated with Newton's inverse square force law, but rather 
with the relative accelerations of neighboring particles (their geodesic deviation), that is, 
with tidal forces. We might say that spacetime curvature and tidal forces represent the 
same physical entity, or to put it in Bondi's words,-- "Einstein's theory of gravitation has 
its physical and logical roots ... in the existence of Newtonian tidal forces." 

From this point on the Newtonian tidal tensor will be referred to as the tidal tensor, 
implying, if not stated otherwise, the classical one. 

To become better acquainted with the procedure of the next sections it is be useful to 
recover the main concepts associated with the most typical tidal perturbation - that induced 
by the spherical potential of a point mass^ $p™(r) = —GM/r, where M is the mass of 
the point mass and G the gravitational constant. By applying Eq. ([3]) the tidal tensor 
in Cartesian coordinates for an arbitrary frame of reference does not look very elegant, 
because once we fix the coordinate origin, the tidal tensor is generally non-diagonal for 
any test particle in a gravitational force field. Because it is a real symmetric matrix, there 
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must be a frame of reference for which r^™ is diagonal.— If we choose a frame of reference 
where the x-axis passes through the particle's position, corresponding to the radial spherical 



r 
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(6) 



From Eq. (jlj) we obtain 



2GM 



dx 



d^ty 'yy 



rZdy 



dFC = r^dz 



GM 



GM 



dy 



dz, 



(7a) 
(7b) 
(7c) 



which represents the cubic dependence of the radial force, and the orthogonal longitudinal 
and latitudinal differential tidal forces, respectively. 

The radial tidal force represents a change in the absolute value of the force vector, and 
the longitudinal and latitudinal vectors represent an angular variation of the components of 
the central force: 

2GM (9F(r) 
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yy 



pm 

zz 



dr 
GM 



Fir) 



(8a) 
(8b) 



The sign has an important meaning. The positive sign (in this case that of the radial tide) 
represents a tensile force and the negative one represents a compressive effect.— 
From Eq. fITal) the radial tidal force induced by a point mass is to first order in x 

2GM 



3-Ax, 



(9) 



where Ax is the distance of the particle undergoing that tidal force from the center of mass 
of the perturbed body, and rem is the distance of this center of mass from the tide raising 
point mass of mass M. Equation ([9]) is a common textbook expression for tidal forces. From 
Eq. ([9]) we recover the celebrated Roche's limit in the usual way 



Axr, 



oche 



Mr^ \ 1/3 



2M 



(10) 



The quantity AxRoche, also called the tidal radius, represents the distance at which the 
disruptive radial tidal force on a system of free particles (a perfectly fluid body) of mass 
Mem exactly opposes the system's self-gravity, F(r) = — GMcm/Ax^, the force that holds it 
together. Inside Roche's limit this body will be torn apart by tidal forces. 

Historically the description of tidal fields beginning with tensor-like objects was less 
successful than its potential or force field counterpart because a tensor field needs more 
entities to be visualized (for example, at each point imagine the vectors representing the 
columns of a matrix), while the concept of a scalar or a vector field is intuitively easier to 
grasp. But the price paid might be worse because textbooks sometimes present lengthy and 
complicated calculations to recover Eq. ([7j). 



A. The tidal field of the spherical heterogeneous mass distribution 

As a good example of tidal fields inside matter distributions, we study the potential of 
a heterogeneous spherical mass distribution. We can show^ that the resulting spherical 
potential is given by 



•-max 



$«ph(r) = -47rG'[- /" p(r')(r')W+ [ " p{rydr' , (11) 

where p(r) is the radial mass density and -Rmax is the maximum extension of the spherical 
mass distribution (if we assume that no cutoff radius exists we can set -Rmax = oo). Simple 
geometric considerations tell us that the radial dependence of the enclosed mass in a sphere 
of radius r and mass density distribution p(r) is 

M{r)=ATT p{r'){r'fdr\ (12) 
Jo 

and therefore we write Eq. (fTT]) as 

$^ph(r) = - - AtcG / p{r')r'dr'. (13) 



r 

The first integral term on the right-hand side of Eq. (ITT]) is the gravitational potential, 
whose derivative is in accordance with Newton's second theorem, which states that the 
gravitational force on a body placed outside a spherical shell of mass M is equivalent to the 
force produced by a point mass at the center of that shell.™ The second integral accounts 
for the effect of the external mass distribution. We do not necessarily need to know the 



radial mass dependence M(r) to calculate r. The potential $ or the mass density function 
p(r) alone can give us all the required information. 

The application of Eq. to Eq. (ITT]) is not so straightforward because the former has 
been defined in Cartesian coordinates and the latter in spherical coordinates. We expect 
that the calculations become easier if we rewrite the tidal tensor in spherical coordinates for 
a spherically symmetric potential. We can transform the tidal tensor, which is the negative 
Hessian operator on a scalar function, into the same negative Hessian matrix in spherical 
coordinates. (This kind of transformation requires the use of tensor algebra, covariant 
differentiation, and a somewhat complicated and long calculation with differential geometric 
methods as shown in the appendix.) The tidal tensor for a spherical potential is 



^sph 



/_^!^ ^ 



r or 

V 

\ r or J 



(14) 



In spherical coordinates the trace elements of t'^p*^ highlight the distinction between a tidal 
force induced by a radial change of the magnitude of the force field and that due to an 
angular one 

g2^sph(^) _ gF(r) 
V - W^~~d^ ^^^^^ 
sph sph 1 d^^^\r) F(r) 

The first and second derivatives with respect to r of Eq. fflTl) are 

d^'^\r) AttG r 2 , GM{t) 

=^J^ )(r ) rfr = — ^ (16) 

^2$sph^^) r , ^ / ^ 2GM{r) , ^ , . 

= ^ / p{r'){r'f dr' + ^^^Gp{r) = + AnGpir). (17) 

or r'^ Jq r'^ 

Newton's second theorem for the force field of a spherical potential emerges automatically 

5$sph(^) GM{r) 



dr r 



2 



(18) 



The direct application of Newton's second theorem is sufficient to determine the force field, 
and no explicit evaluation of the potential derivatives is necessary. However, when dealing 
with tidal forces in the presence of local matter where ^^^^ 7^ 0, we cannot avoid the explicit 
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evaluation of both derivatives and we obtain 
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(19) 



Unlike gravitational force fields it is important to keep in mind that tidal fields depend 
on the inside and the outside matter distributions. If we define pir) = as the enclosed 
medium matter density, that is, the homogenous density that would result from a mass M(r) 
enclosed in the sphere of radius r, we can also write 



^sph 



4nG 



\ 



p(r) - p(r) 















1-, 

'3 



(20) 



Thus we can also say that the strength of the radial tidal forces is proportional to the 
deviation of the local density from the enclosed medium matter density p(r). 

As expected, the trace of r^P'^ is the scalar —4TTGp{r), according to Poisson's equation. 
However, the tidal tensor furnishes more information than Poisson's equation alone, because 
it also gives the individual contributions along the three coordinate axes. The components 
of r'^P'^ depend on the local mass density. The angular variation in d6 or d(f) leads to a corre- 
sponding directional change of the tidal force vector, which, in a spherical mass distribution, 
is a movement along iso-density lines. 

III. THE PHENOMENON OF THE TIDAL FORCE REVERSAL: THE SPHERI- 
CAL MASS DISTRIBUTION 



A. The tidal field for a general mass distribution model and for some special cases 

The spherical heterogeneous mass distribution shows an interesting possibility that cannot 
occur in the traditional case of a point mass or for tidal fields outside mass distributions. 
For special conditions the tidal radial component T^f^ in Eqs. ( IT4|) . (fT9l) . or (l20l) . can change 
its sign. Moving radially inside the mass distribution, the tensile radial tidal force (the 
force responsible for the usual tidal bulges) can become weaker, eventually become zero, 
and even change its sign, that is, become a compressive radial tidal force. For the spherical 
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mass distribution this possibility occurs when the local density nears, equals, or becomes 
greater than the enclosed medium matter density. In the next sections we will show how 
this behavior is a feature of tidal forces inside mass distributions for many geometries. 

The gravitational force never changes sign, but its derivative can, and that can make a 
difference for tidal forces. It all depends on how the local matter density function behaves: 
in some physical conditions we may well have an increase of gravitational forces moving 
away from the central regions. 

We now study this aspect more closely for the spherical heterogeneous mass distribu- 
tion. In particular, let us examine the sign reversal of the tensor component t^^"^ = 
— 9^$'^P'^(r)/9r^ = dF(^r)/dr, and see what it represents. From Eq. fll7p we see that a 
transition from a tensile to a compressive radial force occurs when 

r pr 

- / Pir')ir'f dr' - ATiGp{r) = 0. (21) 
Jo 

An immediate solution is the case of the singular mass density function, p(r) oc 1/r, 
which leads to the exact equality of Eq. (12T|) for any radial distance r. Equation (12T1) implies 
that a spherical mass distribution with a mass density function proportional to the inverse 
of the distance from its center exactly cancels the radial tidal forces everywhere. For any 
radial movement the increase (decrease) of the mass enclosed in a sphere of radius r is such 
that it determines an increase (decrease) of the force field that exactly balances its decrease 
(increase) because of this radial change: the net gravitational force field is constant inside the 
whole mass distribution. The potential $(r) oc $o + cr, with $o = "^(0) and c an arbitrary 
constant, the force F = d^{r)/dr oc c = constant, and there is no gradient of the force field, 
d'^^/dr^ = 0, that is, no tidal forces anywhere. In principle, we could imagine multiple star 
systems, Oort clouds, or entire star clusters moving inside such a mass distribution without 
being affected by a tidal perturbation. 

The mass density function p(r) oc 1/r is interesting as a theoretical limit but is unrealistic 
or at least not very common in astrophysical contexts.— It is a theoretical limiting case 
between a situation that gives rise to tensile or compressive tidal forces. 

For a mass distribution function p(r) oc l/r"^, Eq. (12T]) becomes 

which is satisfied if and only if 7 < 1. Inside a spherical mass distribution with 7 < 1, the 
radial tidal forces are always compressive. The usual tensile radial tides with their induced 
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tidal bulges to which we are accustomed no longer exist. One example is the homogeneous 
mass distribution for 7 = 0, p(r) = po = constant, for which we obtain 

g2,|,sph(^) 4^ 

^rr = = -^Gpo. (23) 

In this case it is not the force field that is constant (it grows linearly with r), but the radial 
tidal component, which is a compression force. 

For 7 < 1, M(r) oc r" with a > 2, that is, the rate of increase of the inner mass 
enclosed in a sphere of growing radius r is faster than the inverse of the square root law. If 
gravitational forces increase with r, the derivative has the opposite sign from what we are 
accustomed to seeing for external or pointlike masses: it is this fact that leads to the tidal 
reversal which is nonexistent for any field external to a mass distribution. We will see that 
it is a common occurrence in the central regions of globular clusters, elliptic galaxies, and 
galaxy bulges. Moreover, the notion of Roche's tidal radius is meaningless. Moving closer 
to this tidal reversal zone, the tidal radius diverges toward infinity, and beyond this limit 
it does not exist. This analytic divergence of the tidal radius should be taken into account 
when we consider stellar mass distributions, which if neglected can lead to numerical errors. 

B. The tidal field for the isothermal mass distribution model 

So far we have analyzed some special cases that are unlikely in the real world (though not 
completely unphysical). We will now consider cases with 7 > 0,^^ where the central density 
does not go to infinity. 

An interesting case of a real spherical heterogeneous mass distribution is that of globular 
clusters.— The surface brightness of these objects is well approximated by King's light curve 
law.— Therefore, if we assume that the visible bright matter represents the real matter 
content (globular star clusters have a very low dark matter content if any). King's surface 
brightness can be recovered by a three-dimensional isothermal mass distribution^ 

Piso(r) = T-fV^ = T^, (24) 
1 + [r jrc) 1 + ?7"^ 

where po is the central density, a core radius (defined as the radius where the surface 
brightness is 50% of the central luminosity and is typically about the order of a parsec in 
globular clusters and from 10^ to 10^ pc for galactic bulges), and 77 = r/rc is a scale length.— 
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The distribution piso is of order with 7 ~ in the central regions and 7 ~ 2 for the 
external ones. Thus, from what we have seen in Sec. IIII At we can expect a tidal reversal 
phenomenon close to the central region. 

We substitute Eq. fl2^ into Eqs. fITT]) and f[T2]) and obtain the potential and the enclosed 
mass at the radius r for the isothermal mass distribution. If we scale all formulas by the 
scale length 77 = r/rc, we obtain 



/ \ ^ ^ 2 r^-arctanr] 1 ( I + rf \ 
^Un) = -^rcGrlpo ^ '- - - log ( t—^ ) ■ (25) 



We use Eqs. (fT2!) . (fT^ . and (fT8|) for the gravitational force field and the tidal tensor elements 
and obtain 

^ (77 — arctan??) , , 

FUn) = -^TTGporM (26) 



A-^rn 2(r/-arctanr/) 1 . 
T„ [7]) = AttGpo (27) 



rl^Tiv) = rt;{v) = -4vrGpo^^^^^^. (28) 

We see that lim^_>o ^isoiv) = —2-nGrlpQ \og{l+r(^^-^ is finite if we assume there is a cutoff 
radius -Rmax = '"c'7max fo^ the dimension of the globular cluster. Similarly, Miso(^7max) = 
47rpo'"c(^?max — arctani^max) is finite and represents the cluster's total mass. There is no 
singularity for Eq. (!26|) because lim^^o(^ ^ arctanr^)/?]^ = 0, which is consistent with the 
fact that gravitational forces must disappear at the center of a spherical mass distribution. 
Note how for large 77 it converges toward an inverse force law instead of an inverse square 
one. 

What do the gravitational and the tidal force fields look like in an isothermal mass 
distribution? In Fig. 1 we show plots of Eqs. ( 1271) (continuous line) and ( l28l) (dashed line), 
both normalized to AnGpo, versus the scale length i]. The gravitational force field increases 
with distance from the center until ~ l.Sr^. The radial tidal field is correspondingly negative, 
which means that there are radial forces of compression in this internal region. From Fig. 1 
we see that there is a relatively intense negative peak for the radial tidal field in the central 
part compared with the positive maxima at about ~ 2.7//. From this behavior it is clear how 
in the central regions of an isothermal sphere the radial tidal compression is much greater 
than the tensile ones. From Eqs. (127|) and (!28|) we find that all the tidal tensor elements at 
the center of the cluster (that is, for 1] 0) have a common value of —AnGpo/S. 
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Compared with the positive maxima at about rj ~ 2.7, which does not exceed ~ 0.03 x 
AnGpo, it is clear how the compressive radial tidal effects in a globular cluster's center might 
be more than ten times greater than any radial tensile tidal force anywhere else. Moreover, 
it is clear that there must be a small region inside such a mass distribution at about ~ 1.5 
scale lengths where the radial tidal effects are zero, beyond which the negative tides take 
over abruptly. 



C. The tidal field for the Plummer sphere mass distribution model 

We briefly repeat the same calculations for the Plummer sphere modeU^ which is some- 
times used as a first approximation for the description of the internal structure of galactic 
bulges. Plummer's model is interesting because it describes the simplest spherical potential 
that does not diverge in the central regions as 

$Pi(r) = - (29) 

where Mt is the total mass. To recover the mass density function, ppi(r), we need to apply 
Poisson's equation directly 



v^»p.M = ;3|:(''^^) = *''Gpp.W- (30) 



Hence 



/^Pl(^) = (l/^°2)5/2 ' (31) 

with po = 3Mj'/47rr|?. By using Eqs. ( IT2l) . (fT^ . and (ITSl) we see again that the gravitational 
force field and the tidal tensor elements are 

^Pl(^) = -^^^/^O^^ (1+^2)3/2 (32) 

4 2/7^ — 1 

r^A^) = 3^^Po ^^^ ^2)5/2 (33) 

r^e\v) = rZiv) = "^^^0^^:^^. (34) 

Plots of the tidal force fields are shown in Fig. 2. Its similarity with the isothermal mass 
distribution is evident, but the gravitational force field is peaked more toward its center and 
consequently there is a much more marked abrupt tidal reversal. Therefore it is expected 
that in galactic centers the effect is less smooth than in globular clusters. 

13 



IV. CONCLUSION AND POSSIBLE DIRECTIONS OF FUTURE RESEARCH 



We have seen that expressing tidal forces in terms of a tensor formahsm naturally leads 
to understanding some features of a tidal field inside mass distributions. As an example 
of its possible uses, we applied it to some astrophysical relevant cases. We saw how an 
inverse radial tidal compressive effect emerges in the central regions of an isothermal model 
and a Plummer model, which are first approximations for globular star clusters and galactic 
bulges mass distributions.— In these models the intensity of the central compressive tidal 
effects depend only on the central density po, and the dimension of this central region is 
of the order of the core radius. The most intense compressive tidal effects are expected 
for globular clusters or galactic bulges with a relatively big core radius and a high central 
density. 

These compressive tides may lead to dynamical effects (as might have occurred in the 
evolutionary history of stellar clusters) in the form of gravitational tidal shocks,^ in tidal 
heating,- in star forming clouds,^ in binary or multiple star systems, in the Oort clouds of 
comets,— and is also expected in open star clusters when their orbital path went through this 
internal region. For example, Oort clouds are dynamic systems that can be very sensitive 
to the external galactic tides.- Galactic radial compressive tides can reduce the Oort cloud 
volume considerably if it is inside the radial galactic tidal reversal zone. But this compression 
would not lead to any gravitational collapse because the Oort cloud comets can be considered 
as a non-self-interacting system of particles. Very different and perhaps dramatic effects 
would cause a volume reduction of an interstellar cloud leading to a gravitational collapse 
and to corresponding star forming process. There is now sufficient evidence that suggests 
how some still to be determined, yet very intense star forming processes occur in the center 
of galaxies up to the present. The tidal reversal phenomenon that occurs in the core radius 
regions might be part of a possible explanation, despite the existence of massive black holes 
whose disruptive tides would prevail only a few light years from its center. (Therefore, 
we would expect a ring shaped star forming region.) Also, the transient and abrupt tidal 
compressions we have seen at about the core radius might be considered as a source for a 
Jeans instability.— 

Isothermal models are also believed to approximate the dark matter halos of galaxies. 
A galactic dark halo has a much lower central density (the order of ~ 0.1Mq/pc^) than 
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those found in globular clusters or galaxy centers (~ 1000Mq/pc^), and therefore any tidal 
contribution of a halo can be neglected. It is tempting to conjecture that during the first 
cosmological phases, where the non-baryonic dark matter density in the universe was much 
higher, if not dominant over bright matter, a compressive instead of a disruptive tidal effect 
could have played a role. This effect might have ignited a star formation process in some 
regions of the universe much sooner than previously expected, as has been observed at high 
redshifts. These lines of research need further study. 
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APPENDIX: THE HESSIAN OPERATOR FOR A SPHERICAL COORDINATE 
SYSTEM 

The nature of the Christoffel symbols necessary to calculate covariant derivatives can be 
found in many textbooks. However these textbooks rarely cover the case of the generalized 
fully covariant connection coefficients that are necessary to obtain the Hessian in spherical 
coordinates. This need is an opportunity to recapitulate covariant differentiation in a broader 
context than usually given in introductory courses on differential geometry. 

If {xi,X2,X3} represent the Cartesian coordinate system and 
{Mi(a;i, X2, X3), M2(xi, X2, X3), M3(xi, X2, xs)} represent a new one, then we might ex- 
pect that by applying the usual transformation law for a second rank tensor, the Hessian 
matrix of a scalar function <l>, Hap = d'^^/dxadxjs, in these new set of coordinates should 
become^^ 

H'^ - _^^sdx,dxs_ ^^^^^ 

duadup dua dup 

Unfortunately, in general d'^^/duadup is not a tensor because this transformation does not 

take into account the fact that the Hessian is the second derivative of the gradient, which is 

not a vector but a covector. As it is known from tensor analysis, to differentiate a covector 

requires covariant differentiation. Thus, to recover the tidal tensor (the negative Hessian of 

the potential) for other coordinate systems, we need to first calculate the covariant derivative 
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of the gradient (written also in the new coordinates). 

The covariant derivative (usually labeled by the symbol ";") of a covector T is 



Tp-.y — Tp^y - F^^^Tq,, (A. 2) 

where T^^^ is the normal "comma derivative" on the component 7 of T, and T"^^ are the 
connection coefficients of the second kind. For Cartesian coordinates all the r°^^ vanish and 
the covariant derivative becomes the usual one. In this case the covector is the gradient of 
$ in the new coordinates = d^/du-y. 

To calculate the connection F^^ we need to distinguish between coordinate bases (also 
called holonomic bases) and non-coordinate bases (anholonomic bases). It is beyond the 
scope of this paper to give a detailed and rigorous description of the differential geometric 
meaning of these complications, and the interested reader is read the relevant literature (for 



example, Ref. Il4l. Sees. 8.4 and 8.5). However, we can clarify the difference briefly as follows. 

Take the mathematician's viewpoint that tangent vectors and directional derivatives are 
the same thing, u = d^- Let u and v be two vector flelds and define their commutator as 

[u, v] = [du, 5v] = 5u(9v - <9v5u. (A.3) 

For any basis {^a} we define the commutation coefficients Cajs'^ by 

[e„, e^] = Cap^ey. (A.4) 

By definition a basis is a "coordinate basis" or holonomic if all the c^^"^ = and a "nonco- 
ordinate basis" (anholonomic) if at least one c^^"' 7^ 0. 

For noncoordinate bases we calculate the connection coefficients as F"^^ = g^^T ^p^, with 
F^^^ the fully covariant connection coefficient of the second kind given by 

where cp^a = 9af,Cp^^- 

For Cartesian coordinates 63, = ^, 6^, = ^, 6^ = ^, all commutation coefficients vanish, 
that is, it is a coordinate basis, and Eq. (lA.SP simplifies to the usual Christoffel symbols. 
For spherical coordinates 

d _1 d _ 1 ^ 
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Ref. y obtain (p. 207) 



—c 



ef 



— c- 



r 

cote 
r 



(A.7) 
(A.8) 
(A.9) 



The connection coefficients for the spherical noncoordinate basis given by Eq. (lA.SP are 

1 



f-e 



f(p 



-r 



r 
1 

r 

cot 6* 



(A.lOa) 
(A. 10b) 
(A.lOc) 



and all other coefficients vanish. 

To calculate the second derivatives of $ (the Hessian), we must calculate the covariant 



derivatives ( 1A.2I) . If we write if/3^($) 
to derive, then 



, with T = V$ the covariant vector 



with the coefficients ( lA.lOl) in the basis (1A.6I) . We obtain 



52$ 



1 1 9$ 

r dOdr 09 

1 1 9$ 



i/g,-($) = (V<l>), ^ - r"..(v$) 



r</> J, gjj^ ^ (90c}r sin 6 

1 1 9$ 



+ -- 



92$ 



cos 9 9$ 



sin 6* dcpdO sin^ 6^ 
1 ia$ cot 9$ 



+ 



+ 



sin^ (902 r dr 90 
The Hessian matrix in spherical coordinates is written as 



(A.ll) 

(A.12a) 
(A. 12b) 
(A. 12c) 
(A.12d) 
(A.12e) 
(A.12f) 



i/($(r,0,0)) 



9f* 
8r2 



r.2 flfl 



1 



92$ 



-r2 ^fl 



1 

r dddr 



y.2 QQ2 ~ r dr 



1 



1 



1 a* 



cos 9 9<I> 



y r sin drd4> 



1 a* 



1 c>2# 



cos 6* 9$ 



sin y Ctp 
92$ _|_ 1 a$ _|_ cot 6> 8$ 



(A.13) 
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Note that the trace of H equals the Laplacian of $ in spherical coordinates. 

In general, Eq. flA.131) is not a particularly useful or elegant expression for the Hessian. 
However, if there is spherical symmetry, this matrix simplifies and the tidal tensor, which 
is the negative Hessian of the gravitational potential $, r = —H, applied to a spherical 
potential ^^^^ becomes the diagonal matrix (1141) . as expected (all angular derivatives ^ and 
^ in Eq. (lAl3l) vanish). 

If we are interested in describing tidal fields in elliptic or disk-like galaxies, we can repeat 
the above calculation to obtain the Hessian in elliptic and cylindrical coordinates respec- 
tively. 
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(1984). 

C. Porciani, A. Dekel, and Y. Hoffman, "Testing tidal-torque theory: I. Spin amplitude and 
direction - Testing tidal-torque Theory: II. Alignment of inertia and shear and the characteristics 
of proto-haloes," MNRAS 332, 325-351 (2002). 

In some of the literature the tidal tensor is identified only with the pure shear deformation, that 
is, the traceless part of Eq. ([3]), Tij — ^Tij, and Eq. ^ is called the "deformation tensor." In 
general, tidal deformations can be a combination of both shear and transverse stretch. We will 
maintain the convention where both stretches are considered because, as we will see, it is more 
consistent with the tidal tensor of general relativity. 

When considering a generic, Cartesian, or radial coordinate system we will write ^(R), 
<&(x,y,z), or <I>(r) accordingly. 

J. Binney and S. Tremaine, Galactic Dynamics (Princeton University Press, Princeton, 1988). 
C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (W. H. Freeman, San Francisco, 
1973). 

H. Bondi, "The foundations of the theory of gravitation," Acta Physica Polonica B 30 (10), 
2859-2862 (1999). 

We will use forces and potentials per unit mass. 

From a purely algebraic point of view, we should in principle calculate the eigenvalues of r 
because every n x n real symmetric matrix is orthogonally similar to a diagonal matrix whose 
diagonal elements are its eigenvalues. Fortunately the diagonalization, at least in this case, is 
much simpler. 

Some authors prefer to use the opposite sign. It is a matter of convention. We follow this 
convention because it emerges directly from the calculations without ambiguity. 
R. Caimmi and L. Secco, "Around the Newton's theorem for ellipsoidal two-component figures," 
Astr. Gesellschaft Jenam Abstract Series, 18, 211 (2001). 
'^^ Newton tried, without success, to show that this theorem also holds for homeoids (ellipsoidal 
shells of matter), but the rigorous proof came about two centuries later, Interesting applica- 
tions and generalizations have been done recently also in Ref. [19 . 

On the other side lim^^o p{f) = cc does not make this case unphysical if we think of the 
existence of a singularity as a black hole at the center of the star cluster, and especially notice 
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how the mass does not diverge, because limr^o M{r) = lirrir^o (27rcr^) = 0. 
The case 7 < on an interval in r impUes an increasing density with r which, even if not 
physically impossible (for example ring-shaped matter distributions), is unlikely, or at least 
much less frequently in nature, and we will not address it here. 

A similar case has been studied by Valluri,'^ but without transforming the tidal tensor to spher- 
ical coordinates. We will compare this profile with Plummer's in Sec. IIIICl 
This density distribution corresponds to a stable arrangement of a self-gravitating system that 
satisfies the equations of hydrostatic equilibrium of an ideal gas at a constant temperature. 
I. R. King, "Surface photometry of elliptical galaxies," Astrophys. J. 222 Part 1, 1-13 (1978). 
M. Masi, "Dynamical effects of the radial galactic tide on an Oort cloud of comets for stars 
with different masses and varying distances from the galactic center," |astro-ph /0403599 [ This 
preprint does not cover the case of the inverse tidal effects. 

The Jeans instability occurs when the thermal pressure of a cloud of stellar gas, which causes 
the cloud to expand, is no longer large enough to balance the gravitational pressure that causes 
the cloud to collapse. In the scenarios we describe here we must consider not only the self- 
gravitational pressure, but the tidal pressure must also be taken into account. 
The models we have considered here are only approximations; more realistic mass distribution 
might give more precise answers. However, the use of these approximations does not undermine 
the fact that the tidal reversal is not an exception for a specific model, but is a universal feature 
of the gravitational field. 
Einstein's summation convention is used. 
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FIGURE CAPTIONS 
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FIG. 1: The tidal force field of the isothermal spherical mass distribution T^°{r]) (continuous line) 
and T^^°(?7) = Tf?{ri) (dashed line). 
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FIG. 2: The tidal force field of the Plummer spherical mass distribution r^^iv) (continuous line) 
and T^^{ri) = t^Uij) (dashed line). 
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